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Abstract. The open clusters' fundamental physical parameters are important tools to under- 
stand the formation and evolution of the Galactic disk and as grounding tests for star formation 
and evolution models. However only a small fraction of the known open clusters in the Milky 
Way has precise determination of distance, reddening, age, metallicity, radial velocity and proper 
motion. One of the major problems in determining these parameters lies on the difficulty to sep- 
arate cluster members from field stars and to assign membership. We propose a decontamination 
method by employing 2MASS data in the encircling region of the clusters NGC1981, NGC2516, 
NGC6494 and Mil. We present a decontaminated CMD of these objects showing the member- 
ship probabilities and structural parameters as derived from King profile fitting. 
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Introduction 

Despite the large number of known open clusters many of them have been discovered 
recently and have not been studied in detail or do not have their paramet ers determined 



(DAML02 - iDias etal. Il2002l: iMermilliod Ill995l : iKharchenko et al. ifiooHh . Sky surveys 
like 2MASS ([Skrutskie et al. I l2QQ6h produced large amounts of near-IR data and have 
contributed to the discovery of even more objects. 

The determination of the fundamental parameters of open clusters provide impor- 
tant tools to the investigation of the Galactic disk and star formation models. However, 
contamination by field-stars and strong reddening severely hamper the detection and 
characterization of clusters, specially towards the Galactic plane or the Galactic center. 

We devised a method to remove the field-star contamination on colour-magnitude dia- 
grams (CMD) by sampling the population from a nearby control field and then removing 
it from the cluster's CMD. The 2MASS catalog was selected owing to its whole sky 
coverage and freedom to extract data from spatially unlimited regions. Also, near-IR 
wavelengths are particularly sensitive t o discriminate cluster st ars from the contaminat- 
ing field for young stellar systems (e.g., Santos Jr. et al. ll2QQ5h . 



Data 

Vizier (http://vizier.u-strasbg.fr/viz-bin/VizieR) was used to extract near-infrared pho- 



tometric data from 2MASS in circular fields centered on objects NGC 1981, NGC2516, 
NGC 6494 and Mil; with the clusters' center coordinates taken from DAML02. The data 
comply with the 2MASS Level 1 Requirement (namely J < 15.8, H < 15.1 K < 14.3 
mag.) and was extracted within 30 arcmin of the cluster's center coordinates. 

Comparison fields were selected for the targets with the same area of the cluster data 
and similar reddening, as deduced from near-infrared (2MASS) and mid-infrared (IRAS) 
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images. NGC 1981 comparison field was extracted 1 deg. northwest due to the southward 
nebulosity associated with NGC 1977. For the remaining targets circular annular fields 
centered in the cluster were used. Figure [1] compares spatial density profiles for the targets 
built from raw data and the decontaminated subsample, discussed below. 




Figure 1. Spatial density profiles for targets (left to right) NGC 1981, NGC 2516, NGC 6494 and 
Mil using raw data (top) and decontaminated data (down). Color bars on the right represent 
stellar density. 



New Center Determination 

Precise determination of the clusters' central coordinates is an essential step before a 
reasonable radial density profile can be obtained. Catalogued center values are mainly 
intended for identification purposes and generally imprecise. Therefore, we have estimated 
the center of the targets by first selecting a region around the coordinates given by 
DAML02. Furthermore the selected region was divided in bins of right ascension and 
declination and star counts were made inside them. We used these star counts to build 
spatial profiles and fit a gaussian function to the star's distribution on RA and DEC. 
The center coordinates were initially taken as the center of the fitted gaussian. 

This procedure was applied to different bin sizes and the corresponding center coor- 
dinates were used to create histograms showing the most recurrent RA and DEC. No 
trends were found between the coordinates and bin size and we adopted these most re- 
current values as the new center coordinates of the cluster. Table [T] show the results of 
the applied method. 

Table 1. (left) Center determination and their difference to the cluster center in DAML02 
catalog, (right) Determined King parameters and uncertainties. 
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Membership Probability Assignment 

For membership assignment we built CMDs for both cluster and field stars and divided 
the diagrams into small cells in J and J — H axes. The cells are small enough to detect 



S266. Membership and CMD Decontamination 121 




Figure 2. Cluster and control field CMDs. From top left to bottom right: NGC1981; NGC 
2516; NGC 6494; Mil. Cluster raw data (left panel), control field (middle panel), cluster cleaned 
subsample (right panel). Color-coded inset show the assigned membership probability. 



local variations of field-star contamination on the various sequences in the CMD, but 
large enough to accommodate a significant number of stars. Typical cell sizes are A J — 
0.5 mag. and A(J — H) = 0.25 mag., with dense clusters accepting smaller cells while 
the sparse ones requiring larger ones. Membership probabilities are assigned to cluster's 
stars within each cell based on the overdensity of cluster stars with relation to the field 
stars, according to the relation P = (N c i u — Nfid)/N c i u . Null probability was assigned 
whenever an excess of field stars over cluster stars occurred in a given cell. 

For the decontamination of the CMD, a subset of the original cluster sample was 
created by removing from each cell in the cluster's CMD, the expected number of field 
stars as measured in the control field CMD. Stars where removed based on their distance 
to the center of the cluster, and cells without cluster overdensity had all stars inside 
their limits removed. We will refer to this final subsample as clean subsample hereafter. 
Fig. [1] shows the spatial density profiles of the targets using the raw data and the clean 
subsample. Fig. [2] shows the cluster, field and clean CMD for the target clusters. 



Radial Density Profile 

In order to probe the spatial profile and derive structural properties of the targets, 
radial density profiles (RDP) were built by counting stars inside rings of one arcmin 
width, centered in the cluster and then dividing by the area of the ring. King profile 
fitting used the modified function cr(R) = ab g -\-cro/(l-\-(R/Rc) 2 )j as introduced in lKing 
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Figure 3. King profile fit on the target clusters. From top left to bottom right: NGC 1981, 
NGC2516, NGC 6494 and Mil. Core radius and background density are indicated. 

(1962). Fig.[3]shows RDP around the new center coordinates for the target clusters using 
the clean subsample. 

Table [T] shows the central surface density <jo, core radius Rc and background surface 
density a^, with their respective uncertainties, for the decontaminated s amples. The val- 
ues fo und for M 11 agree, within the uncertainties, to those determined bv lSantos Jr. et al. 
( 20051 ). 

Conclusion 

By removing stars similar to the field population in the cluster region, the method 
effectively leaves the final subsample of stars less contaminated by the background. Ad- 
ditionally it provides membership probabilities for all cluster stars, removing stars with 
null probability from the sample. This subsample allows for a much better defined cluster 
sequences in the CMD, providing good conditions for subsequent isochrone fitting and 
more accurate values for cluster's parameters, as can be clearly seen in Fig. [T] and Fig. [21 
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